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We have investigated transport through graphene Andreev interferometers exhibiting reentrance 
of the superconducting proximity effect. We observed a crossover in the Andreev conductance 
oscillations as a function of gate voltage (Vbg)- At high Vbg the energy-dependent oscillation 
amplitude exhibits a scaling predicted for non-interacting electrons, which breaks down at low Vbg- 
The phenomenon is a manifestation of electron-electron interactions, whose main effect is to shorten 
the single-particle phase coherence time r^. These results indicate that graphene provides a useful 
experimental platform to investigate the competition between superconducting proximity effect and 
interactions. 



The superconducting proximity effect (PE) consists in 
the modification of the electronic properties of a normal 
conductor (N) in contact with a superconductor (S), by 
induced pair correlations. It is responsible for a broad 
variety of transport phenomena, such as the the modi- 
fication of the tunneling density of states in the normal 
metal [T , the occurrence of Josephson supercurrents [2] 
and subgap structures [3 in superconducting junctions, 
the enhancement of the subgap conductance in single NS- 
junctions [4|, and the modulation of dissipative transport 
by the superconducting phase in Andreev interferome- 
ters [5] . Many of these phenomena are by now well un- 
derstood in terms of Andreev reflection at the N/S in- 
terface, in conjunction with the coherent propagation of 
non- interacting electron and hole waves in N [6], with 
an overall very good agreement between theory and ex- 
periments. On the contrary, except for specific problems 
(e.g., quantum dots connected to superconductors [7 ), 
transport in the regime where electron-electron interac- 
tions (EEI) compete with the superconducting PE has 
remained largely unexplored. 

Here we investigate phase-modulated transport 
through a diffusive Andreev interferometer whose nor- 
mal region is a T-shaped graphene ribbon. We exploit 
the possibility to electrostatically tune the transport 
regime in the ribbon to investigate the competition 
between the PE and Coulomb interactions. Moving 
from large carrier density towards charge-neutrality we 
observe the transition from a regime consistent with 
expectations based on a non-interacting electron picture, 
to one in which the PE is increasingly suppressed. On 
the basis of the experimental observations, we conclude 
that the suppression is due to EEI, whose main effect 
appears to be the shortening of the single-particle 
phase-coherence time r^ in the N region. Our results 
indicate that graphene-based hybrid devices provide an 
excellent platform to explore the PE in the presence of 
interactions and disorder. 

Our investigations focus on a specific manifestation of 
the superconducting PE, namely its so-called reentrance 
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FIG. 1. (Color online) (a) False color scanning electron mi- 
crograph of one of our devices; graphene and the supercon- 
ducting electrodes are colored in cyan and red, respectively 
(contacts 1 and 2 are joined by a :^ 12 pm^ superconduct- 
ing loop), (b) Two-terminal conductance between contacts 
3 and 4 as a function of Vbg, for different values of perpen- 
dicualr magnetic field, (c) Square resistance of the T-shaped 
graphene ribbon versus Vbg- (d) Oscillations in the four- 
probe conductance Gs,2\i,4^ ^^^ ^^ quantum interference of 
Andreev-reflected holes, modulated by the flux threading the 
loop (data measured at T = 250 mK with an applied bias 
VsD =40pV) . 



[8 . This counter-intuitive phenomenon consists in the 
non-monotonic energy dependence of the conductance of 
a diffusive normal metal connected to a superconducting 
electrode through a highly transparent contact. When 
the temperature T is lowered from just above the criti- 
cal temperature Tc, the conductance first increases, and 
then unexpectedly decreases so that at T = it returns to 



the normal-state value Gat, as if the PE was completely 
absent. A similar non-monotonic trend is also observed 
at low temperature, when measuring the differential re- 
sistance as a function of applied bias V (i.e., decreasing 
V from V > A/e to 0; A is the superconducting gap), 
and in Andreev interferometers, when looking at the am- 
plitude of the conductance oscillations as a function of 
bias or temperature. 

Theory treating the normal conductor in the diffusive 
limit relates the energy dependence of the conductance 
change SG{E) to the Thouless energy Et = hD/L'^ 
(where L is the length of the N region, and Et <C A) 
and to the normal state resistance Rn = ^/Gn [S]- 
For a fixed device geometry, the hallmark of the non- 
interacting theory is the prediction of a universal scal- 
ing of the phenomenon in terms of reduced variables, i.e. 
when RnSG{E) is plotted as a function of E/Et 19|. Pi- 
oneering experiments have demonstrated the reentrance 
of the PE in systems where the normal conductor was 
either a thin metal film [lOl [11] or a two-dimensional 
electron gas (2DEG) hosted in an InAs or InGaAs-based 
heterostructure [121 US] • Despite quantitative deviations 
between theory and experiments (especially in the case 
of InAs-based systems) attributed to non-ideal aspects 
of the devices [14 , the excellent qualitative agreement of 
the experiments with theoretical predictions has led to 
the conclusion that the non-interacting theory accounts 
well for the key aspects of the phenomenon. However, the 
universality of the scaling between RnSG{E) and E/Et 
has never been verified experimentally. 

As compared to devices used in the past, graphene 
Andreev interferometers offer two main advantages to ad- 
dress the role of interactions. First, the stability and gate 
tunability of graphene allows us to compare the depen- 
dence of the reentrance effect on the transport parame- 
ters of the normal conductor. Second, in graphene rib- 
bons EEI, enhanced by disorder, become relevant as the 
Fermi energy approaches the charge neutrality point [15 : 
scattering at the edges increases the tendency of elec- 
trons towards (Anderson) localization, which enhances 
the effect of Coulomb interactions (indeed, in sufficiently 
narrow ribbons, fully developed Coulomb blockade is ob- 
served [15 ). It is this latter point that provides a handle 
to tune the effect of EEI experimentally. 

Figure 1(a) shows a SEM micrograph of one of our 
Andreev interferometers. A single-layer graphene flake 
is patterned into a T-shaped ribbon, connected to a su- 
perconducting loop, to control the relative phase of the 
superconducting order parameter Sep = 27r$/<I>o {^ is 
the magnetic flux threading the loop and ^o = h/2e 
the superconducting flux quantum). The device is fabri- 
cated on graphene exfoliated onto a degenerately doped 
Si wafer (coated with a 285 nm thick Si02 layer) acting 
as a gate electrode. Two additional probes are placed on 
the wider graphene region below the vertical arm of the 
"T" . The superconducting loop and these electrodes, are 
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FIG. 2. (Color online) (a) progressive averages of the mag- 
netoconductance measured around Vbg = 501/, for Vsd = 
40/iy, as a function of the number N of used to calculate 
the averages (curves are offset for clarity), (b) Ensemble- 
averaged (thick line) versus individual dl/dV curves (thin 
lines) measured around Vbg = 601/, at T = 250mK. (c) Am- 
plitude of the conductance oscillations measured in the range 
200G < B < 500G as a function of the number N of mea- 
surements used to average the data, (d) Thouless energy Et 
extracted from the VeG-dependent resistance of the device. 



defined first (by electron beam lithography, evaporation, 
and lift-off) and consist of a trilayer of Ti/V/Au (layer 
thicknesses are 5/17/5 nm; the critical temperature Tc 
and superconducting gap A are 3.5 K and 530 peV, re- 
spectively). In a second step graphene is etched in an O2 
plasma, through a PMMA mask, to define the T-shaped 
contour. We have realized and investigated several sim- 
ilar interferometers exhibiting analogous behavior, and 
here we present data measured on one of these devices. 

Experiments were performed in a filtered ^He system, 
down to 250 mK. The two-terminal conductance between 
probes 3 and 4, measured as a function of Vbg in the 
presence of a perpendicular magnetic field B is shown in 
Fig. 1(b): the observation of clear half-integer quantum 
Hall plateaus confirms that the device is fabricated on 
monolayer graphene. The resistance across the ribbon, 
measured at 5 = as a function of Vbg (Fig- 1(c)), 
rises by two-orders of magnitude near the charge neu- 
trality point, due to the mentioned enhanced tendency 
towards localization. Finally, Fig. 1(d) shows the con- 
ductance oscillations induced by a small magnetic field 
that modulates the superconducting phase, originating 
from quantum interference of holes Andreev reflected at 



the two different superconducting contacts. 

As the transport properties of the device are deter- 
mined by quantum interference, we need to distinguish 
between contributions of ensemble-averaged (EA) and 
sample- specific (SS) nature in the measured quantities. 
In particular, the SS component of the magnetoconduc- 
tance oscillations [12 can not be neglected at small en- 
ergies, where the EA component is suppressed by the 
reentrance of the PE. Therefore, for each quantity of in- 
terest we averaged measurements for 25 different values 
of Vbg^ stepping the gate voltage just enough to cause 
a change in Fermi energy larger than the correlation en- 
ergy tt'^Et [H]. A plot of Et = hD/L'^ for positive Vbg 
is shown in Fig. 2(d). The diffusion constant D is esti- 
mated from the zero-bias conductivity, assuming a linear 
dispersion for graphene [17] (we take L = 1pm, corre- 
sponding to the distance between the bottom of the "T" 
and the superconducting contacts; note how D decreases 
as the Fermi level approaches charge neutrality, consis- 
tently with the expected tendency towards localization). 

Figures 2(a,b) show the results of the ensemble aver- 
aging performed around Vbg = 60 V, for the magnetic 
field dependence of the linear conductance and for the 
bias-dependent differential conductance. Around B = 
the amplitude of the conductance oscillations measured 
at VsD = 40peV is not much affected by the averag- 
ing process, because at this bias the EA contribution is 
larger than the SS one. On the contrary, at higher mag- 
netic field {B > 15 mT) where the EA contribution is 
suppressed by the broken time reversal symmetry (due to 
the magnetic fiux threading the graphene "T" ) , averaging 
over N traces suppresses the amplitude proportionally 
to 7V-^/2 (see Fig. 2(c)). Similarly, individual dl/dV 
curves are asymmetric and exhibit random bias depen- 
dent features, whereas the EA curve is symmetric (thick 
versus thin lines in Fig. 2(b)). Note how, with decreasing 
the bias from V > A/e, the EA differential conductance 
dl/dV increases in a smooth featureless way, reaches a 
maximum, and then decreases, so that the value at F = 
is comparable to that measured at V > A/e. This is a 
manifestation of the reentrance of the PE. 

Having established the averaging procedure, we now 
look in detail at the EA phase-modulated oscillations at 
large charge density {Vbg = 60 V). The low-field conduc- 
tance oscillations are plotted in Fig. 3(a) at T = 250 mK 
and for Vsd between and 0.55 meV. Similar measure- 
ments have been done as a function of temperature, for 
Vsd = 0. The Vsd- and T-dependence of the peak-to- 
peak amplitude of the first and second harmonic (whose 
presence is particularly clear at low bias) are shown in 
Fig. 3(b) and (c). The first harmonic exhibits clear re- 
entrance both in the bias and in the temperature depen- 
dence, with the oscillation amplitude having a maximum 
at an energy (i.e., either bias or temperature) comparable 
to the Thouless energy {Et ~ 45 peV at Vbg = 60 V, see 
Fig. 2(d)). The second harmonic, on the contrary, shows 
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FIG. 3. (a) Ensemble- averaged Andreev conductance oscil- 
lations measured at Vbg = 60 V and T = 250 mK, for Vsd 
varying from (bottom curve) to 0.5 mV (top; curves offset 
for clarity). (b),(c) Bias and temperature dependence of the 
amplitude of the first and second harmonics of the oscillations 
shown in (a). 



no reentrance. This is expected, because the trajecto- 
ries causing conductance oscillations with twice the fre- 
quency have to Andreev-refiect at both superconducting 
electrodes, and are therefore longer (by approximately 
twice the distance between the S contacts). As a result, 
the effective Thouless energy associated to these trajec- 
tories is significantly smaller (by approximately a factor 
of 5-6) than Et, so that the energy at which reentrance 
would occur for the second harmonic is smaller than the 
lowest temperature reached in the experiment [18]. 

We now analyze the evolution of the energy depen- 
dence of the oscillations as a function of Vbg- Eig. 4(a) 
shows the bias dependence of the EA oscillation ampli- 
tude (first harmonic) for seven different values of Vbg 
between 60 V and 12.5 V. Upon lowering Vbg^ the maxi- 
mum oscillation amplitude decreases, qualitatively in line 
with expectations based on a non-interacting theory, be- 
cause Rn increases (see Fig. 1(c)). The value of Vsd cor- 
responding to the maximum oscillation amplitude, how- 
ever, remains unexpectedly unchanged. Within a non- 
interacting picture, this finding is inconsistent with the 
value of the Thouless energy Et, which changes from 
~ 45peV to ~ lOpeV as Vbg is lowered (Fig. 2(d)). For 
a more quantitative analysis, we look at the data in terms 
of normalized quantities, i.e. we plot RnSG as a function 
of cVsd/Et (see Fig. 4(b)). At large densities, specif- 
ically for Vbg = 60, 50 and 40 V, the rescaled curves 
fall on top of each other, as expected in the absence 
of EEL Starting from Vbg = 30 V, however, deviations 
from perfect scaling become progressively larger as Vbg 




FIG. 4. (a) Bias dependence of the amplitude of the en- 
semble averaged Andreev conductance oscillations (first har- 
monic) measured at T = 250 mK near zero magnetic field 
(— lOG < B < 10 G), for seven different values of Vbg be- 
tween 60 and 12.5 V (as indicated in the legend; charge neu- 
trality is at Vbg = 10.5 V). (b) Same data as in (a), plotted 
in dimensionless units: the curves measured for Vbg = 60, 50, 
and 40 V exhibit a perfect scaling, which breaks down start- 
ing from Vbg = 30 V. Inset: calculated conductance of a 
diffusive NS junction for increasingly shorter values of Lcp/ L 
(oo, black curve, 0.6, 0.4, 0.3, 0.2, orange curve), using the 
linearized Usadel equations, as in TOi". Note how the maxi- 
mum decreases in amplitude and shifts to larger energies with 
decreasing L^. 



is lowered. Specifically, the maximum relative oscillation 
amplitude decreases, and shifts to larger cVsd/Et ratios. 
Both trends observed in the range 12.5 V < Vbg < 30 V 
are in conflict with what is expected for a non-interacting 
system, and indicate the occurrence of a crossover to a 
different transport regime. 

The evolution of the bias-dependent differential con- 
ductance dl /dViVsD) with Vbg (see Fig. 5(a-c)) pro- 
vides additional information. At Vbg = 60 V (top 
panel) Andreev reflection results in a clear conductance 
increase. At higher temperature (3.5 K) the conductance 
enhancement extends to all subgap voltages while at 
low T a conductance dip appears at low bias (i.e. the 
phenomenology of the reentrance effect). Upon lower- 
ing Vbg to 20 V (Fig. 5(b)), the conductance enhance- 
ment at subgap voltage becomes significantly less pro- 
nounced; eventually, for Vbg sufficiently close to charge 
neutrality {Vbg = 12.5 V, Fig. 5(c)) no enhancement 
of dl /dViVsD) is observed, and only a suppression per- 
sists, which occurs on an energy scale larger than the 



superconducting gap. This suppression is what is typ- 
ically seen in low-dimensional systems where dynamical 
Coulomb blockade becomes relevant [19] , and is expected 
in wide graphene ribbons at temperatures for which a full 
transport gap due to disorder-induced Coulomb blockade 
is not fully developed. We conclude that the manifesta- 
tion of EEI in the dl /dV{VsD) curves occurs over the 
same Vbg range for which the maximum amplitude of 
the conductance oscillations decreases rapidly and shifts 
to higher E/Et ratios (see Fig. 4(b)). For these rea- 
sons, we attribute the observed deviations in the scaling 
between RnSG{E) and E/Et to EEI. 

It is known theoretically that Coulomb interactions 
tend to suppress pairing correlations and compete with 
the superconducting proximity effect. In a context very 
closely related to the one discussed here - the so-called re- 
flectionless tunnelling [20], i.e. the manifestation of the 
proximity effect in devices analogous to ours, but with 
the N and S regions that are tunnel coupled - the sup- 
pression of proximity effect by EEI has been analyzed 
theoretically [21]. It was found that the influence of 
EEI on the electronic phase coherence plays an impor- 
tant role. Indeed, the damping and shift towards higher 
energy of the maximum in the Andreev oscillations ob- 
served in our experiments can be accounted for in terms 
of a progressive shortening of the single-particle coher- 
ence time r^. Specifically, a finite value of the phase- 
beaking length L^ = ^/DT^ introduces a cutoff for the 
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FIG. 5. (a)-(c) Ensemble- averaged differential conductance 
curves between 265 mK (blue) and 3.5 K (red), for different 
values of Vbg (60,20, and 12,5 V, respectively). With ap- 
proaching charge neutrality, the conductance enhancement 
due to Andreev reflection visible in (a) is suppressed (b), and 
eventually completely disappears (c). 



penetration of the pair amplitude in N: when L^ drops 
below L, £^0 = hD/L\ takes the role of Et in deter- 
mining the energy-scale of the reentrance. Having a new 
energy scale (next to Et) that becomes relevant explains 
the deviations from scaling on the energy axis. Dephasing 
obviously also explains why the amplitude of the prox- 
imity effect decreases, since trajectories whose length is 
larger than L^ cannot contribute to phase coherent ef- 
fects. The effect of shortening r^ is illustrated in the 
inset of Fig. 4(b). The plots show the conductance of 
a single NS junction, calculated by solving the linearized 
Usadel equations [10] for five different values of the phase 
coherence time: the trends observed reproduce qualita- 
tively the behavior of the Andreev conductance oscilla- 
tions. At the conceptual level our observations appear to 
be in line and to support the conclusions of theoretical 
studies [21 . 

In conclusion, we have used graphene-based Andreev 
interferometers to investigate the competition between 
the superconducting PE and electron-electron interac- 
tions, which manifests itself through a crossover in the be- 
havior of the Andreev conductance oscillations as a func- 
tion of gate voltage. Interaction effects become more pro- 
nounced as the Fermi level is shifted closer to the charge 
neutrality point, and can be interpreted in terms of a 
progressive suppression of phase coherence in graphene. 
Earlier experiments had already shown that graphene 
is a versatile experimental platform for the investiga- 
tion of proximity-induced superconductivity in conduc- 
tors with tunable transport properties [22]. Our results 
now demonstrate that hybrid devices based on graphene 
nanoribbons are also particularly suitable to study the 
interplay between electron-electron interactions and the 
superconducting proximity effect. 
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